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Abstract 

Helium atoms and hydrogen molecules can be strongly bound inside interstitial channels within 
bundles of carbon nanotubes. An exploration of the low energy and low temperature properties 
of ^He atoms is presented here. Recent study of the analogous ^He system has shown that the 
effect of heterogeneity is to yield a density of states N{E) that is qualitatively different from the 
one-dimensional (ID) form of N{E) that would occur for an ideal set of identical channels. In 
particular, the functional form of N{E) is that of a 4D gas near the very lowest energies and 
a 2D gas at somewhat higher energies. Similar behavior is found here for "^He. The resulting 
thermodynamic behavior of this fermi system is computed, yielding an anomalous form of the heat 
capacity and its dependence on coverage. 
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I. INTRODUCTION 



Extensive efforts have been made to study the essentially one-dimensional (ID) behavior 
of gases exposed to the very linear, confining geometries provided by ensembles of car- 
bon nanotubes and other quasi- ID porous materials, such as zeolites and MCM-41 
Thermodynamic and structural measurements are traditional ways to explore such reduced 
dimensional behavior. One of the more direct signatures of the ID character of the adsor- 
bate should appear in the heat capacity C(T) in the temperature (T) regime where phonon 
excitations characterize the dynamics of the adsorbate at high density ^,0]. In that case, 
the power d of T entering the low T relation C(T) ~ T'^ is a function of the effective di- 

n 

mensionality of the adsorbate: d = 1 for linear confinement, = 2 in a monolayer, etc.[5|] 
Another manifestation of d involves the form of C(T) in the classical Boltzmann regime of 
very low coverage and high T; there, the behavior for a system of N particles is 
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In some cases, however, such seemingly straightforward descriptions of the effective di- 
mensionality can be questioned, leading to more careful investigations, with surprises to be 
found. Such was the case in our recent study (henceforth called I) of particles moving in 
a system of interstitial channels (ICs) between three adjacent carbon nanotubes j^UOl. 
This system had previously been studied extensively with idealized models of the adsorp- 
tion geometry, which predict ID properties of the system. B, y, fl y, y, 0, H Q, Q In 

the new study, in contrast, it was found that the density of states for single particles (^He 
or H2) is drastically altered by heterogeneity; some heterogeneity is always present within 
experimental samples of nanotube bundles We found in I that the density of states 

N{E) near the very lowest energy has the form: 



NiE) = a{E - EminME - (E ~ E™„,) (2) 

Here Emm is the lowest energy level present in the set of ICs, 9(x) is the Heaviside unit step 
function and a is a constant that depends on both the experimental sample and the details 
of the dispersion relation of the particles. 

This linear dependence of N{E) on {E — Emm) is qualitatively distinct from the 
{E — Emin)~^^'^ dependence found for a conventional ID gas. Eq. (2) arises (as discussed 

2 



below) when the distribution of transverse energies associated with the heterogeneity is 
convoluted with the ID density of states describing motion parallel to the tubes' axes. The 
behavior of N{E) at very high energy is given instead by the relation expected for a uniform, 
homogeneous ensemble of Njc identical ICs: 
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{E > Ernin) (3) 
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Here, L is the length of each nanotube (assumed constant) and m is the particle mass. Eq. 
3 describes a strictly ID spectrum of particles. In paper 1 we found that between the two 
energy regimes described by Eqs. (2) and (3) there exists another regime in which the density 
of states is nearly constant: N{E) (E — E^nin)^- Thus, the power A characterizing the 
energy dependence N{E) [E — Emin)^ varies from A=l to A=0 to \—-l/2 as E increases. 
Recall that the relationship upon which Eq. 1 is based presumes that the density of states 
is that of a d-dimensional ideal gas, for which A = d/2 — 1. Hence, the behavior of N{E) 
found in I for the nonuniform IC distribution implies that the effective dimensionality of the 
gas is d=4 at the lowest energy, d=2 at intermediate energy and d=l at high energy. This 
variation leads to a corresponding T-dependent power law behavior of C (T) as T increases. 

In I, we examined the consequences of this unusual spectrum for the very low T behavior of 
an ideal bose gas. The calculated properties include a Bose- Einstein condensation (EEC), 
characterized by 4D thermal properties at both very low T and close to the transition 
temperature, evolving into a 2D gas at somewhat higher T and, eventually, to a ID gas 
at very high T. The present paper is directed instead at calculating the properties of an 
ideal fermi gas in this environment, focusing on the low T properties; for this system, of 
course, there is no transition analogous to EEC. However, the behavior of the fermi system 
is sensitive to the forms of N{E) in the ultra-low T and low T regimes. These regimes 
exhibit 4D and 2D behavior, respectively, for the function C{T). 

The next section evaluates the properties of a fermi gas of '^He atoms exposed to this nan- 
otube bundle. Qualitatively similar behavior is expected to occur for a gas of HD molecules, 
but we have not undertaken the corresponding calculations in that case. The mathemati- 
cal description of N{E) in both paper I and in Section 2 assumes that the ICs within the 
experimental sample include some having the lowest possible energies, near £'^i„. Section 3 
considers an alternative situation, corresponding to the case when the particular sample of 
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tubes and ICs does not include those with such low energies; this alternative case leads to 
quite different predictions for C{T) from those discussed in Section 2. Section 4 summarizes 
our results and draws some conclusions about this interesting system. 



II. PROPERTIES OF THE ^HE GAS 

The methodology employed in I is used here to compute N[E) for individual "^He atoms 
(which we assume to interact only with the sorbent). The single particle states are calculated 
by assuming that the atoms move freely in the z direction (parallel to the tubes' axes); in 
this case, the energy is the sum of a transverse energy and a longitudinal energy. The 
transverse energy is derived by numerical solution of the Schrodinger equation for each 
specified combination of radii surrounding an IC. For a given IC, we define a vector R = 
{Ri, R2, R3) with components equal to the values of these neighboring radii. A particular 
experimental sample is represented as a cloud of points in R space. The function /(R) 
is defined as the density of points in this abstract space. Note that /(R) is an extensive 
variable, which integrates over R to yield the number of ICs present in the experimental 
sample: 



The calculations in this paper employ the same model potential used in I for ^He, so the 
only difference in N{E) arises from the smaller mass of the ^He. For single wall nanotubes 
of typical diameter 1.4 to 2 nm, the IC provides a strongly attractive potential that restricts 
the probability density to close proximity to the center of the IC. Because the atoms arc so 
highly localized, the excited states of the transverse spectrum lie far above the ground state 
energy, which we denote ^^^(R); hence, these excited states may be neglected in computing 
the properties of the system that are relevant to low temperature. The total energy E(R,p) 
of a molecule in a particular state is then a function of four "quantum numbers": the three 
radii and the momentum (p) parallel to the IC direction. 



The density of states at energy E is derived by adding up contributions from all of the ICs 
present in a given experimental sample: 
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(4) 




(5) 



N{E) = ^ / dRf{R)6[E - E{R,p)] 



(6) 
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Here, the sum involves the usual quasicontinuous distribution of momenta. By integrating in 
the usual way over this variable, one obtains an expression involving the transverse density 
of states, g{Et): 



In a hypothetical uniform case, all of the tubes have a common radius (i.e. R = 



{Ro,Ro,Ro) = Ro ) so that /(R) = NicSiR - Ro) ; then g{E) = Nic5\E - ^^(Ro)] 



and the result of evaluating Eq. 6 coincides with the ID solution, Eq. 3, with Et{Ro) 
replacing E^nm- In practice, the singular threshold behavior of Eq. 3 is exquisitely sensitive 
to heterogeneity. As a consequence, the low energy behavior of N{E) is dramatically altered 
from this ID form. We omit the details of the calculation, which are discussed in paper I. 
Fig. 1 presents the density of states for ^He in a distribution of ICs identical to that used 
in I for ^He. There is no evidence of any singular behavior because the inverse square root 
singularity in Eq. 3 is removed by the heterogeneity; in this case, N{E) at very low E is 
seen in Fig. 1 to exhibit the same linear dependence (Eq. 1) on energy above threshold as 
was found in I, although the numerical value of the coefficient a differs, as do the values of 
Emin (-446. 79K for ^He) and Rmin (8.71 A ). This similarity is expected because the only 
difference between this case and that explored in I is the mass. The 4D-like, linear depen- 
dence of N{E) arises from convoluting the ID longitudinal spectrum with a function g{E) 
that is 3D-like, derived from those ICs for which R lies within a 3D annular region in R 
space near near R^m (corresponding to the energy domain [E, E + dE]); see I for a detailed 
discussion of this point. 

The thermal properties of the gas are computed (as in I) by ignoring interparticle in- 
teractions. This assumption lacks rigorous justification but it may be appropriate in the 
nanotube environment due to electrodynamic screening of the interatomic interaction jlfil |. 
Also not considered is elastic screening, associated with the tubes' deformation, which has 
not been studied, to the best of our knowledge. In the analogous case of adsorption on 
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a flat surface, tliese effects are usually small 
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18l Il9j . but the present situation differs 



because the adsorbate wave function is localized closer to the nanotubes than it is to the 
planar surface. Indeed, we have found at T=0 that the interaction leads to a non-negligible 



dilation of the lattice o: 
at low density 
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nanotubes at high particle density and a factor of 10 smaller dilation 



With the neglect of all such interactions, the energy and specific heat are computed in 
the conventional way for an ideal fermi gas. For example, the energy of the system satisfies 

U = J dEEN{E){exp[(3{E - + (9) 

Here l/jS = ksT and /i is the chemical potential, which is determined from the equation for 
the total number of atoms: 

N = j dEN{E){exp[f3{E - //)] + 1}-^ (10) 

The classical limiting behavior corresponds to the case of low density, when exp[/3{E — fi)] S> 
1. Then, the energy per particle is 

[/ d{lnQ) 

\ ) classical \ / 

where the single particle partition function is 

Q = J dEN{E)exp[-(3E]. (12) 

The Fermi energy corresponding to the low energy limiting form of the density of states, 
Eq. 2, is given by the relation 

Ef = — (13) 

a 

Fig. 2 shows the behavior of Ef over the energy regime extending to higher energy, the 
so-called 2D regime where N{E) is constant. 

The effective dimensionality of the system varies with both N and T since the low T 
behavior probes the Ef dependence of N{Ef), which varies with N. Figure 3 presents re- 
sults obtained for C{T) = {dU/dT)j\f at very low T. The behavior at the lowest density 
corresponds closely to the classical limit (A^ approaching 0), except for T < ImK, when 
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degeneracy effects appear. That limit has the behavior described quahtatively in the intro- 
duction. The T = classical limit is C/lNks) = 2, since the lowest energy density of states 
has the form (A^ proportional to E) of a 4D gas. Between 1 mK and 15 mK the behavior 
is that of a classical gas with the density of states shown in Fig. 1. At higher T (above 
0.02 K), the result is C/lNks) — 1, that of a 2D gas. At much higher T (not shown), the 
classical specific heat becomes that of a ID gas, C/{NkB) = 1/2, since then the effects of 
heterogeneity are minimal. 

This behavior evolves as the number of atoms, A^, increases. At intermediate values of A^, 
the degenerate regime of C ~ T extends to higher T before converging to the classical limit. 
At the highest density appearing in Fig. 3a, there is little evidence of the 4D behavior and 
the result is similar to that of a strictly 2D gas. This is demonstrated in Fig. 2 by an explicit 
comparison with the prediction for the 2D gas. Note that the low T limiting behavior is 
C/^B = ttT, where a = {n'^ /3)k'^N{Ef) , i.e., proportional to the density of states at the 
fermi level. When the density of states is linear in E, a is proportional to A^^^^. The general 
dependence can be undestood from Eqs. 1 and 2, which determine how a depends on A^. 

III. ALTERNATIVE (GAUSSIAN) MODEL CALCULATION 

In the preceding discussion, the sample's distribution of tubes is assumed to include ICs 
in the vicinity of the global energy minimum, Emin- In this section, we consider the extreme 
opposite case, in which the distribution function /(R) is centered about a region of R space 
near a point Rq that is far from Rmm- We assume a gaussian form for this distribution: 



Here a is the distribution's width parameter and F = A'"/c'(27r)^^/^(T~^ is a normalization 
constant, yielding the integral N/c forf (R). The transverse density of states is obtained 
from Eq. 8, which yields the relation [2^ 



Here, the integral is over the equal energy surface S in R space that satisfies E = Et(R). In 
the following, we assume that a is "small", so that a linear expansion of the energy in the 
vicinity of Eq = £'t(Ro) suffices to characterize the transverse energy within the distribution: 



/(R) = Fexp[-(R-Ro)V(2^')] 



(14) 




(15) 
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Et{R)^Eto + VEf{R-Ro) (16) 

Here Efo = i?i(Ro) is the most probable transverse energy within the distribution of ICs 
and VEt is the gradient evaluated at i?o-With this expansion and the gaussian distribution 
function, Eq. 14 can be evaluated analytically, yielding 

As might have been expected, this distribution assumes a gaussian form. The standard 
deviation of the energy distribution a \ VEt \ is large if the cloud of IC points representing 
the sample is spread broadly in R space and/or the energy gradient is large. 

Our analytical result [17] for this gaussian model represents a stark contrast to the result 
g{Et) ~ {Et—EminY^'^ found near the energy threshold in I, a consequence of the distribution 
of ICs being localized around R„wn ; it is the 3/2 power law dependence that yields the d=4 
low energy spectrum discussed in the previous section of this paper. In the present case, 
Eq. 17 may be integrated, using Eq. 7, to obtain the total N{E). The result appears for 
two sets of parameters in Fig.4 

One observes in these curves that the ID character oi N{E) is preserved at high E, but an 
exponential tail of low energy states appears, which will determine the low energy/T behavior 
of the adsorbed gas. Fig 5 exhibits the heat capacity obtained with varying coverages with 
this model. 

The most interesting feature of these results is the extreme sensitivity of the low energy 
and low T behaviors to the distribution of ICs. At the very lowest energy, N{E) is very flat, 
which might be interpreted as a very high effective dimensionality, meaning that C assumes 
a high value in the classical regime; see Eq. 1. The coefficient of the specific heat in the 
linear, degenerate regime is very sensitive to both the number of atoms present and the form 
of N{E), as usual. 



IV. SUMMARY AND CONCLUSIONS 

In this paper we have explored the spectrum and thermal properties of a gas of ^He 
atoms confined within a set of ICs. Two distinct situations were considered. The first study 
assumed that the distribution of ICs in R space is localized near the global energy minimum, 
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as was assumed in I. In this case, the low T (4D) behavior of C{T) is very dependent, as 
seen in Fig. 3. There occurs a maximum in C{T) during the crossover between the regimes 
of degenerate and classical statistics. One can, in principle, have ID, 2D or 4D behavior in 
classical and quantum regimes. The description, in general, is complicated because the low 
T behavior is determined by the number of particles, which determines the density of states 
at the fermi level. In the alternative gaussian model the low T behavior is very different 
from that of the preceding model. The low T behavior is determined by the tail of the 
exponential at low A^, which does not correspond to any particular effective dimensionality. 
In contrast, at high (in both models) the properties are those of a ID system. 

This research has been supported by the National Science Foundation. We are indebted 
to Mercedes Calbi, Carlo Carraro, Tony Clark, Susana Hernandez, Karl Johnson and John 
Saunders for discussions of this work. 
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FIGURE CAPTIONS 

1. The density of states of ^He atoms in the vicinity of the lowest energy Emin 

2. Fermi energy as a function of A^. 

3. a) The specific heat of a gas of -^He atoms at densities N = 10'^ A'^, N = lO^^A"^ 
and N = lO^^A^'^. b) Specific heat of an ideal 2D fermi gas at the same densities as in a). 

4. Densities of states N{E) obtained from the gaussian model of the nanotube distribution 
in R space, with the parameters cr = 0.1 K (full curve) and cr = 0.3 K (dashed curve). 

5. Low T heat capacity obtained with the gaussian model, with cr = 0.1 K and number 
densities N = 10~^A-^ = lO^^A"^ and A^ = lO'^A"^ (from top to bottom). 
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a = 0.01K 
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